Regularization of Synchronized Chaotic Bursts 
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The onset of regular bursts in a group of irregularly bursting neurons with different individual 
properties is one of the most interesting dynamical properties found in neurobiological systems. In 
this paper we show how synchronization among chaotically bursting cells can lead to the onset of 
regular bursting. In order to clearly present the mechanism behind such regularization we model the 
individual dynamics of each cell with a simple two-dimensional map that produces chaotic bursting 
behavior similar to biological neurons. 



PACS number(s): 05.45. +b, 87.22.-q 

Studies of cooperative behavior in coupled chaotic os- 
cillators are frequently based upon the analysis of the 
phenomenon of chaos synchronization M . Different types 
of synchrony between chaotic oscillators along with the 
various mechanisms responsible for the onset of such syn- 
chronization have been studied 0. These studies are 
mostly motivated by the development of a theoretical 
framework that can explain synchronization in neurobio- 
logical systems where complex, chaotic behavior of neu- 
rons is quite typical 

It has been observed in neurobiological experiments 
and in numerical simulations that individual neurons 
may show irregular bursts, while ensembles of such ir- 
regularly bursting neurons can synchronize and produce 
regular, rhythmical bursting. This regularization occurs 
despite the significant differences in the individual dy- 
namics of coupled neurons [0] . How can a group of dif- 
ferent neurons, whose individual dynamics are very sen- 
sitive to the intrinsic parameters of the neuron, synchro- 
nize and produce a robust rhythm of the bursts? How 
can such rhythmical bursting become insensitive to the 
intrinsic parameters that control the individual dynam- 
ics of the neurons? Answers to these questions may be 
found through the study of models, which, on one hand, 
are capable of demonstrating similar effects and, on the 
other hand, are simple enough to reveal the mechanisms 
that lead to these effects. 

In this communication we present an example that 
demonstrates the effects of mutual synchronization and 
chaos regularization of bursts in a group of chaotically 
bursting cells. We study N oscillators, which are modeled 
with two-dimensional maps and coupled to each other 
through the mean field H 



x(i,n + l) = — — 

(1 + x{i,nY) 

y(i, n + l) = y(i, n) - ax(i, n) - (3 , 



N 



(1) 



where x(i,n) and y(i,n) are, respectively, the fast and 
slow dynamical variables of the i-th oscillator, e is the 
strength of global coupling. The slow evolution of y(i, n) 



is due to the small values of the positive parameters /3 and 
a, which are on the order of 0.001. The value of parame- 
ter a is selected in the region a > 4.0, where the map pro- 
duces chaotic oscillations in x(i, n). When e = 0, depend- 
ing on the value of parameter a, each cell demonstrates 
two qualitatively different regimes of chaotic behavior: 
continuous chaotic oscillations and chaotic bursts, see 
Fig. 0a. Such dynamics resemble the bursting dynam- 
ics measured in the experiments with biological neurons, 
see for example B. 
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FIG. 1. The waveforms (left) and power spectra (right) 
of the chaotic behavior of uncoupled cells, e = 0, computed 
for different values of a, with a = (3 — 0.001. The power 
spectra are computed for waveforms of length L, (L=50000 
iterations) . 

As one can see from Fig. 0a, the mean duration of the 
bursts is very sensitive to the value of a. The chaotic 
nature of the duration of the bursts can be clearly seen 
from the power spectra calculated from x(i, n). The low- 
frequency part of the spectrum that corresponds to the 
frequency range of the bursts is presented in Fig. 0b. 
The introduction of coupling between the cells results in 
the synchronization of the bursts. While the bursts do 
become synchronized, the fast chaotic oscillations of the 
cells during the bursts remain asynchronous. The most 
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interesting effect observed in this case is that the syn- 
chronized bursts measured both in each chaotic cell and 
in the mean field (X}j=i X (j> n )/N) become almost peri- 
odic, see Fig |2^,. We define the formation of such a peri- 
odic temporal pattern in the chaotic behavior of systems 
as chaos regularization. The effect of chaos regulariza- 
tion can be clearly seen in the power spectrum of x(i, n) 
as the appearance of periodic spectral components (com- 
pare Fig. [l]b and Fig. ||b). These periodic components 
correspond to the period and the harmonics of the peri- 
odic temporal pattern formed in the synchronized chaotic 
bursts. In the rest of the paper we will study the mech- 
anisms behind the chaos regularization of the bursts in 
the model cells M). 



c„ + i = F(x n ,a,j) = 
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FIG. 2. The waveforms (left) and power spectra (right) of 
the synchronized chaotic bursts of the coupled cells (hf) with 
e — 0.2 and N = 256. The values of the individual parameters 
of the cells are randomly selected from the intervals: a from 
4.1 to 4.9; (3 and a from 0.0009 to 0.0011. The value of L is 
the same as in Fig. [j] 

Chaotic Bursting. The individual behavior of the cells 
is described by the map of the form 



(1+4) 
Vn+1 = Vn - C.T„ - P 



(2) 
(3) 



The slow evolution of y n for the next m steps is given by 

% (4) 



Vn+r 



where x n ,m — Ej=n+i x i)l m ^ s the mean value of x n 
computed for m consecutive iterations. It is clear from 
(|]) that the value of y n slowly increases during the next 
m steps if <jx n ,m < — 0, and decreases if ax niTn > —ft. 

Since y n changes slowly, the dynamics of x n can be 
considered independently of the map (|^) assuming that 
y n is a parameter 7 = y n . Therefore, fast dynamics of 
the cell can be understood from the analysis of the 1-D 
map 
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The main dynamical properties of this map are clear 
from Fig ||. The results of the bifurcation analysis of 
the map are presented in Fig. ||. If the parameter val- 
ues of the map are within the horn formed by the bi- 
furcation curves L T \2 and L T 23, given by the equations 
a = -2(( 7 2 + 9)7 - ( 7 2 - 3) 3/2 )/27 and a = -2(( 7 2 + 
9)7 + (7 2 - 3) 3/2 )/27 respectively, then the map (g) has 
three fixed points x*, x\ and x\. The fixed point x\ is 
stable, X2 is unstable and a;^ may change stability, see 
Fig. § 
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FIG. 3. The shape of the function F(x„,a,"/) plotted for 
a — 4.1 and three different values of 7 (-2.65, -2.75 and -2.85) 
(a) and the chaotic trajectory of the map (b). a;*, x% and 
X3 are the fixed points of the map. 7* corresponds to the 
bifurcation value where x\ and x\ merge together. 

When the parameter values of the map cross the curve 
L T i2 from left to right then the fixed points x\ and x% 
merge together and disappear. This bifurcation in (|) 
corresponds to the beginning of the burst in system (2) 
and (^). Indeed, when x n is in the stable fixed point x J, 
then 7 = y n slowly grows because x U} i « x\ < —(3 /a, 
see (Q). When the stable fixed point x\ disappears, the 
trajectory of the fast map (||) goes to the chaotic attrac- 
tor corresponding to the chaotic pulsations during the 
burst, see Fig. ||. 

The end of the chaotic burst is due to the external 
crises of the chaotic attractor in system (||). When the 
trajectory of this map is on the chaotic attractor, then 
x„~i > —P/ a (see Fig. ||), and the value of 7 = y„ de- 
creases until the chaotic trajectory of the fast map ar- 
rives back to the stable fixed point x*. This becomes 
possible only after the parameter of the fast map goes 
through the bifurcation value, which corresponds to the 
case when the trajectory from the maximum of the map 
function x max — F(0,a,~/) maps into the unstable fixed 
point x\. Curve L^, that corresponds to the bifurcation 
values where x max iterates into a fixed point x\ or x*, 
is given by the equation a = —(37 ± \/^ 2 — 8)/2 and is 
shown in Fig. || . Note that the external crises of the 
chaotic attractor is associated only with the fixed point 
X2 ■ This bifurcation corresponds to the low branch of the 
curve Lh that starts from point A12. 
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FIG. 4. The bifurcation diagram of the uncoupled fast map 
(a), and the enlarge portion of the diagram around the point 
A12 (b). The bifurcation curves corresponding to the merging 
of the fixed points x* and x* m (curves L Tjm ) and to the crisis 
of the chaotic attractor (curve L^)- g shows the gap between 
the bifurcation curves. 

Therefore, the duration of the chaotic burst is deter- 
mined by the time interval that is required for the slow 
variable y n to move from the bifurcation curve L T i2 to 
the curve Lh plus the additional time interval that is 
needed for the trajectory to map below the unstable fixed 
point x'2 ■ Both of these intervals fluctuate because of the 
chaotic nature of the trajectory x n , and this is the rea- 
son for the chaotic fluctuations of the duration of the 
bursts. After the chaotic firing terminates the trajectory 
arrives at the stable fixed point x\, and 7 slowly moves 
towards £ r i2, where the next burst starts. One can see 
from Fig that the size of gap g between the bifurca- 
tion curves L T i 2 and is very sensitive to the value of 
parameter a. That explains the high sensitivity of the 
mean duration of the bursts to the internal parameter of 
the cell a. 

If the value of parameter a is selected above the point 
A12, then the fixed points x\ and x 2 appear within the 
chaotic attractor (internal crisis). In this case y n fluctu- 
ates in the vicinity of the bifurcation curve L r i2 keeping 
the cell in the regime of continuous chaotic oscillations, 
see Fig. |l|a. We would like to note that the considered 
mechanism of chaotic bursting is essentially the same as 
in the Hindmarsh-Rose model of biological neuron, where 
the role of parameter a is played by a hyperpolarization 
current, see Ref. Q. 

Synchronization. Coupling between the cells (|l|) influ- 
ences the fast dynamics of each cell (|^) by adding the 
value e^2f = ix(j,n)/N to the parameter 7. When the 
i-th cell approaches the bifurcation values at L T i 2l being 
in the stable fixed point x\, its behavior becomes very 
sensitive to the influence of the other cells. The cells, 
which are crossing this curve and starting to fire, sharply 
increase their values of x(J, n). As a result, the increased 
value of 7 in the i-th cell pushes the cell over the bi- 
furcation value and triggers the burst. The greater the 
number of cells which are involved in the simultaneous 
transition to the firing phase, the larger the triggering 
impact will be that is experienced by the remaining cell. 
Such an avalanche type of synchronous switching takes 



place both at the beginning and at the end of each burst. 
This threshold mechanism of synchronization is similar 
to the one studied in ensembles of integrate-and-fire os- 
cillators ||. 
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FIG. 5. The waveforms of the mean field (solid line) and 
of x(i,n) for the cell with a = 4.5 (dots) computed for the 
group of synchronously bursting cells with e = 0.1 (a) and 
e = 0.2 (b). The bifurcation curves of the fast maps in the 
case of synchronously bursting cells (c). 

Regularization The dynamical mechanism of regular- 
ization can be also understood from the bifurcations of 
the fast map (|^) taking into account the alternation 
of the two phases in the waveform of the mean field 
e 2~2j=i x (j) n )/N for synchronized bursts, see Fig. ||a,b. 
When the synchronously bursting cells are in the non- 
firing phase, the fast variable of each cell is located in 
the stable fixed point x* and the values of y(i, n) and the 
cc-coordinates of the fixed point x*(i,n) slowly increase. 
In this case the evolution of x in each cell can be de- 
scribed by the map ([!]), where the coupling term can be 
approximated by the value ex(i,n). 



x(i, n + 1) = 



(1 + x(i,n) 2 ) 



7 + ex(i, n). 



(6) 



The bifurcation curve L T i2, which defines the threshold 
value of y(i,n) corresponding to the beginning of the 
burst, is now given by the equation 27a(l — e) 2 /2 = 
-( 7 2 + 9(1 - e) 2 ) 7 - ( 7 2 - 3(1 - e) 2 ) 3 / 2 . The location of 
this curve on the parameter plane (7, a) depends upon 
the value of the coupling parameter e and shifts to the 
right from the original position of L T i 2 , see Fig. ||c. 

When all synchronously bursting cells are in the phase 
of chaotic firing their fast chaotic pulsations are un- 
corrected. Therefore, the value of the coupling term 
(e 2^2^=1 X U> n ) /N) during the burst can be approximated 
with the mean value of < x(j, n) >. Due to the averaging, 
the fast chaotic fluctuation in the coupling term disap- 
pears and the value \e < x(J,n) > \ slowly drifts in the 
interval between and 0.08, see Figs || a,b. These val- 
ues are significantly less than the amplitude of the slow 
variable \y„ \, which is of the order of one. Therefore, for 
simplicity, the contribution from the mean field to the 
parameter 7 can be neglected during the firing phase. 
As a result, the approximate position of the bifurcation 
curves for the fast chaotic oscillations of synchronously 
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firing cells (|lj) can considered to be the same as for the 
uncoupled map (Bh. 

Fig. ||c presents the bifurcation curve L T i2 (e > 0) 
that defines the beginning of a burst, and the bifurca- 
tion curves L% (the part below ^412) and L T i 2 (e = 0, 
the part above ^412), which define the end of the burst 
in each synchronously bursting cell. It follows from this 
figure that the states of slow variable y(i,n) correspond- 
ing to the beginning and the end of the bursts become 
separated for all values a. This separation defines the 
period of the periodic pattern formed in the chaotically 
bursting cells and helps to synchronize cells with qualita- 
tively different individual dynamics whose bursting rates 
become insensitive to the parameter a. 

Assuming that variables y(i,n) drift with approxi- 
mately the same rate, the mean period of bursts, Tg, 
should be proportional to the length of the most narrow 
gap g between these two curves. One can see from Fig.||c 
that the size of such a minimal gap g m in increases with 
the strength of the coupling. This explains the change of 
bursts frequency in the waveforms shown in FigJ|a,b. 

The presented analysis shows that chaos regulariza- 
tion in a group of chaotically bursting cells has the fol- 
lowing mechanism. When the large enough fraction of 
the group of cells starts to burst simultaneously due to 
the threshold synchronization, it leads to the formation 
of relatively sharp transitions (ramps) in the mean field 
(coupling term). These ramps create a gap between the 
states of y(i,n) corresponding to the beginning and the 
end of the bursts for the rest of the cells, including the 
continuously firing ones. All cells become bursting cells 
and it improves the synchronization among them. Due 
to the formation of the gap, a large portion of the group 
have about the same mean period which is short enough 
to lead the bursting of the whole ensemble, see Fig.|]c. 
During the chaotic burst the chaotic components of mean 
field are small because of averaging over the whole group 
of the cell and, as the result, these components do not 
influence the individual oscillations of the cell. When 
a large number of firing cells get close to the thresh- 
old, the leading cells trigger an avalanche of transitions. 
This avalanche sharply switches the mean field to the 
non-firing phase and the rest of the cells are forced to 
change their state independently of their chaotic trajec- 
tory. This is how the synchronization makes the duration 
of the bursts insensitive to the chaotic trajectory. It is 
clear that the regularity of the bursts improves when the 
number of cells increases. However, due to the formation 
of the gap, this type of chaos regularization becomes no- 
ticeable even for small groups of cells, including the case 
N = 2. 

Finally, we would like to emphasize that this mecha- 
nism of chaos regularization is due to the dynamical fea- 
tures of each cell at the beginning and at the end of the 
chaotic burst. In the considered map these features are 
characterized by bifurcations at L T \2 and Lh, which are 



similar to a saddle-node bifurcation and the appearance 
of a homoclinic orbit in more realistic models of biolog- 
ical neurons, see for example R|. Since the mechanism 
of chaos regularization relies only upon the formation of 
sharp ramps in the coupling forces, and the averaging 
down of the chaotic fluctuation during the firing phase, 
this mechanism can occur for a number of different types 
of coupling among the cells. 
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